Levy flights in quenched random force fields 



Hans C. Fogedby 

* Institute of Physics and Astronomy, University of Aarhus, DK-8000, Aarhus C, Denmark 

and 

NORDITA, Blegdamsvej 17, DK-2100, Copenhagen 0, Denmark 
(February 1, 2008) 

Levy flights, characterized by the microscopic step index /, are for / < 2 (the case of rare 
events) considered in short range and long range quenched random force fields with arbitrary vector 
character to first loop order in an expansion about the critical dimension 2/ — 2 in the short range 
case and the critical fall-off exponent 2/ — 2 in the long range case. By means of a dynamic 
renormalization group analysis based on the momentum shell integration method, we determine 
flows, fixed point, and the associated scaling properties for the probability distribution and the 
frequency and wave number dependent diffusion coefficient. Unlike the case of ordinary Brownian 
motion in a quenched force field characterized by a single critical dimension or fall-off exponent 
d — 2, two critical dimensions appear in the Levy case. A critical dimension (or fall-off exponent) 
d — f below which the diffusion coefficient exhibits anomalous scaling behavior, i.e, algebraic spatial 
behavior and long time tails, and a critical dimension (or fall-off exponent) d = 2/ — 2 below which 
the force correlations characterized by a non trivial fixed point become relevant. As a general result 
we find in all cases that the dynamic exponent z, characterizing the mean square displacement, locks 
onto the Levy index /, independent of dimension and independent of the presence of weak quenched 
disorder. 
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I. INTRODUCTION 

There is a current interest in the dynamics of fluctuating manifolds in quenched random environments This 
fundamental issue in modern condensed matter physics is encountered in problems as diverse as vortex motion in high 
temperature superconductors, moving interfaces in porous media, and random field magnets and spin glasses. In this 
context the simplest case is that of a random walker in a random environment, corresponding to a zero dimensional 
fluctuating manifold. This problem has been treated extensively in the literature and many results are known. 

In the case of ordinary Brownian motion, characterized by a finite mean square step, in a pure environment without 
disorder, the central limit theorem JH) implies that the statistics of the walk is given by a Gaussian distribution with 
a mean square deviation proportional to the number of steps or, equivalently, the elapsed time, i.e., the mean square 
displacement is 

(r 2 {t)) oc Dt 2 / Z , (1.1) 

where the dynamic exponent z assumes the value z = 2 for Brownian walk; D is the effective diffusion coefficient for 
the process and (• ■ •) denotes an ensemble average. 

There are, however, many interesting processes in nature which are characterized by anomalous diffusion with 
dynamic exponent z ^ 2, owing to the statistical properties of the environments J0JJ. Examples are found in chaotic 
systems ||, turbulence (7]||, flow in fractal geometries and Levy flights Jl0|jllj |; these cases generally lead to 
enhanced diffusion or superdiffusion with z < 2; we note that the ballistic case corresponds to z = 1. The other case 
of subdiffusion or dispersive behavior with z > 2 is encountered in various constrained systems such as doped crystals, 



glasses or fractals [12 16 



Independent of the spatial dimension d, ordinary Brownian motion traces out a manifold of fractal dimension oIf = 2 
p7[ . In the presence of a quenched disordered force field in d dimensions the Brownian walk is unaffected for d > dp, 
i.e., for d larger than the critical dimension d c \ = dp the walk is transparent and the dynamic exponent z locks onto 
the value 2 for the pure case. Below the critical dimension d c \ = 2 the long time characteristics of the walk is changed 
to subdiffusive behavior with z > 2 p^-|2Q|. In d = 1, (r 2 (t)) oc [logi] 4 , independent of the strength of the quenched 
disorder 
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Levy flights constitute an interesting generalization of ordinary Brownian walks. Here the step size is drawn from 
a Levy distribution characterized by the step index / [§,[ll). The Levy distribution has a long range algebraic tail 
corresponding to large but infrequent steps, so-called rare events. This step distribution has the interesting property 
that the central limit theorem does not hold in its usual form. For / > 2 the second moment or mean square deviation 
of the step distribution is finite, the central limit theorem holds, and the dynamic exponent z for the Levy walk locks 
onto 2, corresponding to ordinary diffusive behavior; however, for / < 2 the mean square step deviation diverges, 
the rare large step events prevail determining the long time behavior, and the dynamic exponent z depends on the 
microscopic step index / according to the relationship z — / (/ < 2), indicating anomalous enhanced diffusion, that 
is superdiffusion pl] , p2| , p3[ . The "built in" superdiffusive character of Levy flights has been used to model a variety 
of physical processes such as self diffusion in micelle systems [g4| and transport in heterogeneous rocks [g2| . 

In a recent letter |25[ we considered Levy flights in the presence of a quenched isotropic random force field and 
examined the interplay between the "built in" superdiffusive behavior of the Levy flights and the pinning effect of the 
random environment generally leading to subdiffusive behavior. Generalizing the discussion in |l8|,^0|,^6| we found 
that in the case of enhanced diffusion for / < 2 the dynamic exponent z locks onto /, independent of the presence of 
weak disorder. On the other hand, we could still identify a critical dimension d c i = 2/ — 2, depending on the step 
index / for / < 2. Below d c \ the quenched disorder becomes relevant as indicated by the emergence of a non-trivial 
fixed point in the renormalization group analysis. We also briefly discussed the behavior of the subleading diffusive 
term. Here another critical dimension, d C 2 = f enters in the sense that for dimensions less than d C 2 the correction to 
the ordinary diffusion coefficient diverges. This feature has also been discussed from a heuristic point of view in ref. 
. Note that in the Brownian case for / = 2 the critical dimensions coincide, i.e., d c \ = d c i- 
In the present paper we discuss in more detail Levy flights in quenched random force field with arbitrary range 
and vector character thus corroborating and extending the results obtained in ref. | p5| . The paper is organized in 
the following manner. In Sec. II we introduce and discuss Levy flights in a pure environment. Next we derive the 
Langevin equation for Levy flights in a force field and the associated Fokker Planck equation for the probability 
distribution. In Sec. Ill we establish perturbation theory for the probability distribution in the Fokker Planck 
equation averaged over the quenched random force field and set up the renormalization group analysis. In order to 
extract the scaling properties and control divergent contributions we dilute the degrees of freedom by means of the 
momentum shell integration method. Finally we derive differential renormalization group equations to first loop order 
for the parameters in the model. In Sec. IV we determine flows and fixed point structure for the renormalization 
group equations in some detail in the case of short and long range isotropic force fields and determine the scaling 
properties. Since the discussion for anisotropic force fields is quite analogous we summarize this case in the Appendix. 
We conclude the paper with a summary and a conclusion in Sec V. A range of technical issues are for the record 
deferred to an extensive Appendix. 

II. THE MODEL 

It is our aim to analyze the scaling properties of Levy flights in a quenched random force field. For that purpose 
let us first discuss and summarize the properties of Levy flights in a pure environment. 

A. Levy flights 

We consider a particle or for that matter a gas of non-interacting particles performing independent isotropic random 
motion in d dimensions. For the step size distribution we assume a normalized Levy distribution |5|, |rT| , [i^p8| -|3l|l 

p(v)d d V=^V- 1 - f dvdn and S d = _ . (2.1) 

Here r) is the microscopic step and the algebraic behavior in the large step regime is characterized by the microscopic 
step index /. In order to ensure normalization of the distribution p(fj) we introduce a lower cut-off rjo of the order of 
a microscopic length and assume a positive step index, / > 0. Sd is the full solid angle on the surface of a unit sphere 
in d dimensions. In Fig. 1 we have depicted the Levy distribution P(rj). 

The macroscopic physics ensuing from the Levy distribution for the microscopic elementary step depends entirely on 
the range characteristics of P(r]). For / > 2 the second moment or mean square step, (n 2 ) — J p{ri)rfd d n — jz^Vo i 

exists and a characteristic step size is given by the root mean square deviation (n 2 )^ . For 1 < / < 2 the second 
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moment diverges, however, the mean step (rj) = j p(r/)r/d d rj — -jrzjVo is finite, defining an effective step size. In the 
interval < / < 1 the first moment diverges and even a mean step size is not defined. 

For a Levy flight consisting of ./V steps the end-to-end distance of the random walk is — Y^iLi Vi or ; assuming 
one step per unit time, t = N, r(t) = 5Zi=i Vi > an d it follows generally that the mean square displacement for 
statistically independent step events is given by 



(r 2 (t))=t( V 2 ) 



(2.2) 



Consequently, for general step distributions with a finite second moment or mean square step, as in particular the case 
for Levy distribution for / > 2, the mean square displacement is proportional to the elapsed time. This behavior is 
characteristic of ordinary Brownian motion leading to diffusive behavior for a gas of Brownian walker and corresponds 
to the dynamic exponent z = 2 in Eq. (1.1) @J3T 



On the other hand, for Levy flight with step index / < 2 the second moment (r] 2 ) diverges, Eq.(2.2) is undefined 
and we c annot define a mean square displacement and a dynamic exponent z according to the heuristic definition in 
Eq. (1.1). However, this point will be clarified below, in terms of the probability distribution. 

The probability distribution for the random walker is given by P(r, t) = (S(r — r{t))) and is easily derived by means 
of the method of characteristic functions j|,[l7],[30|] . We obtain, noting that t is discrete, 



P(r,t) = J e ife -[p(fc)f 



d d k 
{2n) d 



(2.3) 



where p(k) is the Fourier transform of the step distribution p(rj), p(k) = J e~ lk " n p(r])d d ri . At large r Eq. (2.3) 
samples the small k region. For small k the behavior of p(k) is controlled by the algebraic power law tail of p(rj). We 
find 



p(k) = 1 - D(kr]QY ~ exp [-D^k^] , 



(2.4) 



where I? is a dimensionless geometrical factor and p, is a scaling index. For / > 2 the second moment (i] 2 ) exists and 
ji locks onto 2; for / < 2, (rj 2 ) diverges and p = /. By insertion we obtain 



P(r,t) = J exp [ik -r-\t\ Drfik*]-^ 



(2.5) 



In Fig. 2 we have shown the dependence of the scaling index p on the step index /. From the structure of Eq. (2.5) 
we infer that P(r, t) has the scaling form 



P(r,t)=|*l" G(r/\t\-) 



(2.6) 



For / > 2, ji = 2 and the scaling function G{x) takes a Gaussian form characteristic of ordinary Brownian motion, 
G(x) ~ exp(const.a; 2 ). This is a consequence of the central limit theorem |5j which here implies a universal behavior 
characterized by the exponent p = 2. For / < 2, p — f and the scaling function G{x) can in general not be 
given explicitly in terms of known functions. For p — 1, the ballistic case r ~ t, we find the Cauchy distribution 
G{x) ~ (1 + const. a; 2 )^^ 1 )/ 2 ) . It is, moreover, easy to show that G( x) — » const, for x — > an d G (x) — > for x — > oo. 

The scaling properties of Levy flights are described by Eq s. (2.5) and (2.6). From Eq. ( |2.6| ) we infer the mean 
square displacement (r 2 (t)} ~ t 2 ^, i.e., according to Eq. (1.1) the dynamic exponent z — p, indicating supcrdiffusive 
behavior. However, as mentioned above this reasoning is not correct for Levy flights. The expression in Eq. (|2.2|) 



is not defined for Levy fl ight s since (rj 2 ) diverges and the scaling regime must be defined with some care if we wish 
to give meaning to Eq. (1.1). One way is to confine the Levy flights to a volume of linear size L; the mean square 



displacement is then given by (r 2 (t))L = J r 2 P(r, t)d d r , where V = L d . Using the scaling form in Eq. (2.6) with 



we obtain, changing variables, the mean square displacement 



r 2 ( t )) L =\t\ 2 ^ 



r 2 G{r)d d r 



(2.7) 



(2.8) 



From Eq. ( |2.6| ) and (2.8) we infer that the characteristic time in the problem is the time ~ it takes the random 
walker to traverse the volume V ~ L d . At long times, i.e., for t 1 ^ ^> L, 
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(r 2 (t)) L = \t\ 



/ r 2 G{r)d d r = const. \t 
Jo 



2/n 



(2.9) 



and comparing with Eq. (1.1 
for r large, we obtain 



we deduce the dynamic exponent z = fi. At short times, t 1 ^ <C L, using G(r 



(r 2 (t)) L = const. \t\L 2 



(2.10) 



in accordance with Eq. (2.2). For Levy flights with /i < 2 (r diverges for t — > oo. In the Brownian case 



/.t = 2 and (r (t)) oc |t| in both cases in accordance with Eq. (2.2). Summarizing, in order to characterize anomalous 



superdiffusion arising from Levy flights by means of the mean square displacement we must conceptually confine the 
flights to a box of size L and define the scaling region for t^> LP. Alte rnatively, we can discuss the scaling properties 
and the identification of z by means of the scaling form in Eq. (2.6) where the argument r controls the range in 
question. 



B. Langevin equation for Levy nights 

It is convenient to discuss Levy flights in terms of a Langevin equation with "power law" noise [32,33]. In an 
arbitrary force field F (r), representing the quenched disordered environment, the equation takes the form 

- r ^- = F(r(t))+r,(t) . (2.11) 



Here r/(i) is the instantly correlated power law white noise with the isotropic distribution given in Eq. (2.1) at a 
particular time instant. 

The microscopic steps r)i wit h dist ribution p(rji) constituting a Levy flight are discrete processes, so properly 
speaking the Langevin equation ( [2.11 ) is the continuum limit of the corresponding difference equation defined for a 



discrete time step A. This limit is singular and it is instructive to discuss the limiting procedure in some detail. 

Limiting our discussion to the force-free case, F = 0, setting r n = r(t n ),rj n = rj(t n ), and t n = nA, the difference 
equation is given by {r n+1 — r n )/A = r) n with solution r n = A J2p=o Vv From the definition of P(r, n) = (S(r — r n )), 
we obtain, using S(r) = J exp (— ik • r)d d k/{2ir) d together with exp (ik • r n ) = rj™ =0 exp(zAfc • tj p ) for P(k,n) = 
J exp(ifc • r)P(r,n)d d r, 

P(k,n) = (exp (iAk- 77))" (2.12) 



and the issue is to define the continuum limit A — > 0, keeping t = nA fixed, of the expression for P(k, n) in Eq. 



For simplicity we first consider the case of a Gaussian distribution for p{rf), i.e., the case of ordinary Brownian motion, 
p(v) — (l/2cr7r 1 / 2 ) d exp(— ?7 2 /4cr 2 ) of width 2a. Evaluating (exp(iAfe • r])) by completing the square in the exponent 
we obtain P(k, n) = exp (— A 2 cr 2 fc 2 n). In the continuum limit Act 2 — > D and we have the Gaussian distribution 

P(k,t) = exp[~Dk 2 t] . (2.13) 

Correspondingly, the white noise condition, i.e., uncorrelated noise, (TJnVm) = o^&nm&ap, becomes in the continuum 
limit 

(V a {t)v^{t')) = D8 ap 8{t - t') . (2.14) 

Consequently, in the Gaussian case we must in the continuum limit A — > scale the width a of the noise distribution to 
infini ty i n order to obta in a finite diffusion coefficient D — a 2 A. In the Levy case the noise distribution p(rf) is given by 
Eq. and from Eq.(|3), (exp(iAk ■ rj)) = l-D(kAr] ) f for / < 2 and kA-q < 1, i.e., P(k,n) = [1- D{kAr]o) f } n . 



Using the relation, (1 — x/n) n — > exp(— x) for n — > 00, we obtain in the continuum limit A — > and t = nA fixed, 

P(k, t) = exp [- J D? ? ^A / - 1 /fc / 1] . (2.15) 

In order to eliminate the discrete time step A and keep the coefficient D fixed we must renormalize the cut off r/o in the 
Levy distribution ( [Q| ) according to t]q A-^ -1 = 1. We notice that the border line case is / = 1. For l</<2r;o^oo 
for A — ► 0, i.e., the cut-off moves out to infinity. For < / < 1, on the other hand, 770 —> for A — > 0. We also notice 
that A?7o — > so that kAr] <C 1 is satisfied for all /. 
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Summarizing, we can without loss of generality discuss Levy flights (and Brownian motion) in terms of a continuous 
Langevin equation provided we scale the underlying discrete noise distributions accordingly. Note , how ever, that these 
renormalizations are not observable, the problem at hand is defined by the Langevin equations (2A1). 

Levy flights in an arb itrary force field are in principle described by the Langevin equation ( 2.11 ) together with the 
distribution in Eq. (2.1) for the noise 77. For a given force field the only random aspect resides in the noise 77 which 
drives the position r of the random walker; the force field F(r) acts as a stati c bac kground. However, for a random 
force field modelling the quenched static environment the Langevin equation (2.11) harbours two different kinds of 
stochasticity and it is convenient to recast the problem in terms of an associated Fokker-Planck equation, thereby 
absorbing the fluctuating Levy noise 77(f). 



C. Fokker-Planck equation for Levy flights 



In the absence of the force field, i.e., F(r) = 0, it follows from Eq. (|2 

dP(k,t) 



that P(k,t) satisfies the equation 



dt 



-Dxk»P(k,t) , 



(2.16) 



where we have absorbed the renormalized cut-off 770 in the Levy coefficient D\ — Dt]q. The form of Eq. (2.16) leads 
us to suggest the following Fokker-Planck equation for a Levy walker in a force field: 



dP(r,t) 
dt 



= -V • (F(r)P(r, f)) + fliVP(r, f) + D 2 V 2 P(r, f) 



(2.17) 



Here the first term on the right hand side of Eq . ([2.17 ) is the usual drift term due to the motion of the walker in 
the force field, the second term arises from Eq. ( 2.16| ) where we have introduced the "fractional gradient operator" 



V M as the Fourier transform of — fc M . 



V^ 1 



-a— ju 



is a spatially nonlocal integral operator reflecting the long 



range character of the Levy steps; for fi = 2 it reduces to the usual Laplace operator describing ordinary diffusion 
pijl . Finally, we have for later purposes included the ordinary diffusion term, L>2V 2 P(r, f), originating from the next 
to leading long ra nge part of the distribution ^(77) and corresponding to the next leading term of order fc 2 in the 
expansion in Eq. (2.4). The derivation of the Fokker-Planck equation is discussed in more detail in Appendix A. 
For the quenched random force field we assume a Gaussian distribution, 



p(F(r)) ocexp[-i J d d rd d r' F a {r)A af) (r - r 
with the spatial correlations given by 

{F a {r)F p {r')) F = A afi {r-r') 



(2.18) 

(2.19) 
f denotes 



Here A a/3 (r — r') is the force correlation function expressing t he ra nge and vector character of F(r) and 
an average over the force field according to the distribution ( 2.1 8| ) . 

In the unconstrained case A Q/3 is diagonal; however, generally the force field breaks up into a longitudinal curl-free 
part E and a transverse divergence- free part B, i.e., F = E + B, where V B = and Vx£ = fl Qj2^,^]. Assuming 
that the cross correlation of E and B vanishes, (E a B 13 ) — 0, we have in Fourier space 



(E a {k)E l3 (k')) F = {2n) d 8{k + fc') 
(B a (k)B f3 (k')) F = {2n) d 5(k + fc') 



k a k? 



k 2 



A L (fc) 

k a k 



k 2 



A T (fc) 



(2.20) 
(2.21) 



The case of an unconstrained force field then corresponds to A L = A T = A and we obtain the force correlation 
function 



(F a (k)F^{k')) F = {2n) d S(k + k')5 aP A{k) . 



(2.22) 



The range of the force correlations is characterized by the functions A L (fc) and A T (fc). 

In the case of a finite range, e.g., A L ' T (r) oc exp(— rfco) jr d ~ 2 , with range parameter 1/fco we have A L ' T (k) oc 
l/(fc 2 + fc 2 ). In the long wavelength limit k — > A i,T (fc) — > const., corresponding to the zero range case A L ' T (r) oc 
5(r). For infinite range the force correlations are characterized by the exponent a, A L ' T (r) oc r~ a , A L - T (k) oc k a ~ d . 
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Summarizing, we characterize a general force field by the following expressions for the longitudinal and transverse 
force correlation functions: 

A L (k) = Af + A%k aL - d (2.23) 
A T (fc) = Aj + A?k aT - d . (2.24) 

We have here for completion introduced two separate indices, and ax, for the power law behavior of the longitudinal 
and transverse correlations, respectively. For a^,T > d, A L,T (k) — ► A^' T , in the long wavelength limit and we 
effectively retrieve the zero range case; for a^ x < d,A L,T (k) — > Aj ' T k aL • T ~ d and the specific long range behavior 
enters in the k — > limit. 



III. RENORMALIZATION GROUP THEORY 



The problem of analyzing the asymptotic long time scaling properties of Levy fligh ts in a random quenched force 
field has now been r educe d to an analysis of the random Fokker-Planck equation (2.17), in conjunction with the force 
distribution in Eq. ( 2.18 ) and the force correlation functions in Eqs. ( [0£ - f2~2l" ). In particular we wish to evaluate 
the scaling behavior of the distribution P(r,t) averaged with respect to the force field, i.e., (P(r,t))F- 



A. Perturbation theory 



There are a variety of techniques available in order to treat the random Fokker-Planck equation (2.17). Applying the 
Martin-Siggia-Rosc formalism in functional form pq-|4C[| and using either the replica method M or an explicit causal 
time dependence |18 3^|h|, one can average over the quenched force field and construct an effective field theory. A 
more dire ct me thod, which we shall adhere to in the present discussion, amounts to an expansion of the Fokker-Planck 
equat ion (2.17) in powers of the force field and an average over products of F(r) according to the distribution in Eq. 

(Hi- 

Defining the Laplace- Fourier transform 



P(k,(j) = J d 3 rdtexp(iujt - ik ■ r)P{r,t)0{t) 



(3.1) 



where 9{t) is the step function, we obtain, introducing for later purposes the dimensionless coupling strengths \l and 
At for the vertices coupling to E and B, respectively, the Fokker-Planck equation 



P(k, w) = G Q (k, w)Po(fe) + G (fc, u) J -lP[(-i)A z k • E(k -p) + (-i)X T k ■ B(k - p]P(p, u) 



(3.2) 



Here the force field is averaged according to Eqs. (2.20- 2.21 ) for all pairwise contractions according to the Wick's 
theorem following from the Gaussian distribution in Eq. (2.18), Po(k) = P(k,t = 0) is the initial distribution, and 
we have, moreover, introduced the unperturbed propagator or Green's function 



G (fc,w) 



-ico ■ 



Dik» + D 2 k 2 



(3.3) 



The integral equation ( p^ ) immediately lends itself to a direct expansion in p owe rs of E, B or Xl, At- In order to 
discuss the various perturbative contributions it is convenient to represent Eq. ( |3.2| ) diagrammatically as done in Fig. 
3 ||j2^,^l|-Q. Here P(k,oj) is characterized by a solid bar, Po(k) by a cross, the vertices and — i\r by dots, 

the force fields E and B by wiggly lines, and G (k,ui) by a directed arrow. 

Hence, iterating Eq. ( |3.2[ ) in powers of Aj, and At and averaging over the force fields, keeping one component 
fixed, we identify perturbative corrections to 1) the self energy, 2) the vertex functions, and 3) the force correlation 
functions. Defining the self energy E(fc,w) by means of the Dyson equation, 



G{k,uj) = G Q (k,u) + G (fe,w)E(fc,w)G(A;,w) 



(3.4) 



shown diagrammatically in Fig. 4, where the renormalized propagator is indicated by a solid directed line and the 
self energy by a circle, we derive the renormalized Fokker-Planck equation shown in Fig. 5, 



G 



P(k, w) = G(k, u)P (k) + G(k, u) J 0^ [(-i)A L (k,p, to) ■ E(k -p) + (-i)A T (k,p, u) ■ B(k - p)} P(p, w) . (3.5) 

Here Ai(fe,p, lo) and AT(fc,p,w), depicted as circles in Fig. 5, are the renormalized vertex functions; to lowest order 
A-L,T{k,p,Lu) — X^ xk. In a similar manner we extract corrections to the force correlation function from the 4-point 
vertex function T(k,p, 1,uj) in Fig. 6 and the contraction j T(k,k,l,Lo)G(k — l,uj)d d l/(2ir) d , also shown in Fig. 6. 
To lowest order T(k,k,l,Lo) = -k a A al3 {l){k - if, summed over a and (3, where A a @ is the force correlation function 
defined in Eqs. ( 2.1 9j - p721 ). To second order in Al and At or first order in A L and A T , corresponding to first loop 



order in the field theoretic formulation |18|J36| -|40(| we find diagrammatic contributions to £, Al,t, and T shown in Fig. 
7, Fig. 8, and Fig. 9, respectively. 



Let us first discuss the self energy correction. Solving the Dyson equation (3T) the self energy £(fc,w) enters in 
the renormalized propagator, 

G(fc ' w) = ^JTw+W^somo ' (3 - 6) 

and directly determines the diffusional character of the random walker. In Appendix B we discuss in some detail the 
evaluation of Y,(k,Lo) to leading order in k 2 on the basis of the diagrams in Fig. 7. 

In the case of an isotropic unconstrained zero range force correlation function, Al = A^ = A and A L (k) = A T (k) — 
A, we have in particular 

S(fe, lo) = —A 2 A j ^k ■ (fc/2 + p)G (k/2 +p,u). (3.7) 

To leading order in k 2 t he static contribution, S(fe, 0), given in Appendix B, only contributes to the ordinary diffusion 
term D 2 k 2 in Eq. ( |3.6| ); there is no correction to the anomalous Levy term D\k^. We shall see later that this has 
a profound effect on the scaling properties of Levy flights. For the correction to the diffusion coefficient D 2 we then 
find, performing the integration over the solid angle, 



SD 2 = l -X 2 A^ [ A dp W~»Y +d - l +W-^ d+l f (3 . 8) 
" 2 d(2n) d J 1 [D lP n + D 2 p 2 } 2 



where Sd is given in Eq. ( ]2 . l| ) and we have introduced a UV cut- off corresponding to a microscopic length of order 
1/A. In the long wavelength limit p — > the integrand in Eq. ( |3.8| ) is dominated by the leading Levy term ~ p M and 
simple power counting shows that the integral is convergent for d > \x yielding a correction to D 2 - For d < fi the 
integral diverges in the infrared limit p — > and we need a renormalization group approach in order to disentangle the 
true asymptotic scaling behavior. We encounter here the first critical dimension d c2 — fi characterizing the behavior 
of the ordinary diffusion coefficient D 2 . In the Brownian case p8|-p0t [i = 2 and the critical dimension is 2. 

In a similar way we can discuss the vertex correction A(fc,p, oj) on the basis of the diagrams in Fig. 8 in the 
isotropic case. The vertex function A describes the coupling of the random walker to the quenched force field. In 
Appendix B we discuss in detail the evaluation of A. We obtain 

A(k,p,u) = Xk- A 3 A J -0^{k - l)k • (p - l)G (k - l,o;)G (p - I, to) (3.9) 
from which we extract to leading order in k and for lo = the perturbative correction to A, 



s d r A P d+1 

'd(2ir) d J P [D lP » + D 2P 2 



Similarly to our discussion of the self energy, the integral in Eq. ( 3.1 C| ) is convergent for d > 2/i — 2, whereas for 



d < 2/i — 2 the correction 5X diverges in the far infrared limit p — > 0. We note here the appearance of a second critical 
dimension d c \ = 2/j, — 2, characterizing the behavior of the vertex correction. In the Brownian case [L = 2 and both 
the vertex and self energy diverges for d < 2. We also note that for // = 2 the first order correction to D 2 actually 
vanishes for d = 2 thus requiring an expansion to second order in A (fourth order in A) or, equivalently, to two-loop 
order P|- 

Finally, we discuss the correction to the force correlation function A extracted from the contraction of the 4-point 
vertex function depicted in Figs. 6 and 9. From the results in Appendix B we deduce 
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T(k,p,l,w) 



k- (p- l)X 2 A 

d d n 



-A 4 A 2 



+A 4 A ; 



(27r) d 
d d n 



k • (p — n)(k — n) ■ (p — l)Go(k — n,uj)Go(p — n,i 



k • (p — l)(k — n) • (p — I + n)Go(k — n, u)Go(p — I + n 7 ui) 



and we obtain, contracting Eq. (3.11), to leading order in k and for lu = the perturbative correction to A, 



SA = A A 



2 \ 2 °d 



(2tt)« 



i r p 



(3.11) 



(3.12) 



Also here we note that the correction to A diverges for d < d c ±. In the Brownian case, /i = 2, the critical dimensions 
coincide, i.e., d c \ = d C 2- 



B. Momentum shell integration 

In order to disentangle the breakdown of primitive perturbation theory and deduce the scaling properties of the 
force averaged distribution (P(r,t)) F and the mean square displacement ((r (t ))) F we carry out a dynamic renor- 
malization group analysis, following the momentum shell integration method |^,^l|,^3|,^5) . This approach is a way of 
systematically diluting the short wavelength degrees of freedom, keeping the long wavelength modes controlling the 
asymptotic scaling behavior. The method derives from Wilson's original momentum space procedure [ phi applied to 
static critical phenomena |42f and is an implementation of the real space Kadanoff construction in momentum space 
p3[ . The application of renormalization group theory to dynamic phenomena described by Langevin-type equations 
was initiated in the context of dynamical critical phenomena in refs. ]4^j48| . The implementation of the momentum 
shell integration method was introduced and discussed in refs. pl]-|43[]. 

Here we briefly review the momentum shell integration method in the context of the Fokker-Planck equation ([3.2]), 
which in the case of an isotropic force field is expressed in the symbolic form 

P = G P + XG FP . (3.13) 

Setting the UV cut-off A = 1, dividing the wave number interval < k < 1 into a long wavelength regime < k < 
e~ e , I > 0, including the long distance scaling region, and a short wavelength "shell" region e~ e < k < 1, the idea 



is no w to average out the short wavelength degrees of freedom, using the distribution of the force field in Eqs. (2.18- 



2.19) in the shell e < k < 1, in order to derive a renormalized Fokker-Planck equation valid for long wavelengths 



< k < e . Projecting the Fokker-Planck equation (3.13) onto the two regions in wave number space we obtain the 
coupled equation 

P' = G' Q P^ + XG' (FP + F'P + FP' + F'P') (3.14) 
P = G P + XGq(FP + F'P + FP' + F'P') , (3.15) 



where the prime refers to the shell e < k < 1. Averaging the "long wavelength" Fokker-Planck equation (3.15) with 
respect to the short wavelength degrees of freedom in the shell we have, noting that (F')f' = 0, 

(P) F ' = G Po + XG FP + XG F(P') F , + XG (F'P') F , . (3.16) 

The evaluation of (P')f' and (F'P') p> to a given order in A is now achieved by expanding the "short wavelength" 



Fokker-Planck equation (3.14) and averaging over F' . Using that GqFG' and GqF'G' F'G' vanish in the long 
wavelength limit, since two small wave numbers cannot add up to a wave number in the shell, and defining the force 
contraction according to the notation F' C .F' C = (F'F')dij, we obtain 

{P) F , = G P + XG F(P) F , + X 2 G G F' C G' F' C (P) F , + X 4 G F^G' F^G' F^G' F^(P) F , 

+X i G F' Ci G' Q F' C2 G' Q F' C2 G' F' Ci (P) F , + X 3 G F^G' FG' F^P) F , + X 4> GqF' c G'qFG' FG' (s F' c (P) F i (3.17) 

which is diagrammatically depicted in Fig. 10. 

We note that the momentum shell integration method combined with a perturbative expansion is essentially a 
nonlinear procedure leading to a more general Fokker-Planck equation involving higher order force fields as shown in 
Fig. 10 where the last diagram corresponds to two force fields F coupling to the distribution P. More importantly for 
deriving renormalization group equations we identify the same self energy, vertex, and force corrections as in primitive 
perturbation theory but now evaluated within the short wavelength shell. In this manner the elimination of the short 
wavelength degrees of freedom enters in the Fokker-Planck equation for the remaining long wavelength degrees of 
freedom. 
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C. Renormalization group equations 



While the general expression to first loop order for the self energy, vertex, and force correlation corrections and 
the general renormalization group equations are given in Appendix B, we here present a detailed derivation of the 
renormalization group equations in the iso tropi c short range case. 

Disregarding the two- force term in Eq. ( |3.17 ) which is of higher order in k and therefore becomes irrelevant in the 
long wavelength limit, we obtain the renormalized Fokker-Planck equation 



£>!*" + (D a + 5D 2 )k 2 ]P(k, w) = P (fe) + (A + SX) J t^4(-*)* ■ F(k - p)P(p, i 



and force correlation function 



(F a (k)F^p)) F = (A + SA)S a ^{2n) d S{k + p) , 



(3.18) 



(3.19) 



both valid for wave numbers in the interval for < k,p < e~ l . 

The first step in deriving renormalization group equations is then to rescale the wave number range to < k,p < 1 
and in this manner compensate for the eliminated degrees of freedom. Rescaling at the same time frequency, probability 
distribution, and force field, according to 

k' = ke l 



I al£) 

uj = cue y ' 



(3.20) 
(3.21) 

P'(k',u') = P{k,u)e- aW (3.22) 
F'(fe') = F{k)e~ m , (3.23) 
where a(£) and [3(£) are to be determined subsequently, we thus obtain the renormalized Fokker-Planck equation 

[-iu' + D[(£)k'^+D 2 (£)k' 2 }P'(k',oj')=P^k') + (-i)X , (£) J ^tL k ' -F'(k' -p')P'{k'^) (3.24) 
and force correlation function 

(F°"(k')FP'(p')) F , = A'{£)S al3 (2Tr) d 5{k' +p') . 



for wave numbers k' and p' in the original range < k',p' < 1; note that Po(k) 
the scale dependent parameters D[(£) , D' 2 (£) , \'(£), and A'(^) given by 

D[(£) = D ie - ifi+a ^ 
D' 2 (£) = (D 2 + SD 2 (£))e~ 2e+a W 
\'{£) = (A + 5\{£))e-^ 1+d ^ a ^ + ^ 
A'(£) = (A + <5A(<'))e d ^ 2/3 W . 



(3.25) 

P(5(fc'). We have here introduced 

(3.26) 
(3.27) 
(3.28) 
(3.29) 



Note that the corrections SD 2 (£),SX(£), and SA(£) depend on the scale parameter £ since they are evaluated in the 
short wavelength shell e~ t < k < 1. 

In order to allow for an iteration the renormalization group equations are re cast in diff erential form by considering 
an infinitesimal scale parameter £ and expanding the right hand sides of Eqs. ( 3.26| -3.29). Defining 



a{£) 



z(e')d£' 



and noting from Eqs. (3^ 



/3(e) = / u{£')d£' 
Jo 

( 3.10) ), and ( 3.1 2| ) evaluated on the shell k,p = 1 that 



6D 2 oc 
SX oc 



1 S d lD 1 (d- ( jL) + D 2 (d-2) 

2 (27r) d d 
S d 1 



[Dx 
A 3 A 



X 2 A£ 



(2Tr) d d[D 1 +D 2 } 2 
A 2 A 



£ 



(2tt) c 



d'[D 1 + D 2 f ■ 



(3.30) 
(3.31) 

(3.32) 
(3.33) 
(3.34) 
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We arrive at the differential renormalization group equations 

d -^ = {z- i x)D 1 (3.35) 

— = (, - 2)D 2 + A [Di + D2? A A (3.36) 

d\ A 3 A 

5-c+- 1 -*-V^ (M7) 

§ = (( ,-2„)A + 2^-l) 15 ^ p , (3.38) 

where A = {l/2d)S d /(2n) d is a geometrical factor, £ d = 2?: d / 2 /{d/2 - 1)! (cf. Eq. Q). 

In conformity with Ref. p(| we have included a vertex coupling A. However, since A 2 is always associated with A 
in the diagrammatic expansion, the inclusion of both A and A is essentially superfluous and we can for examp le se t 
A = 1 and discuss the coupling to the force field by means of A alone Thus assuming d\(£)/d£ = in Eq. (3.37) 
and solving for u(£) we finally obtain the renormalization group equations 

-jf = {z -»)Di (3-39) 

f =( & -< i - 2 ' a - 2 ' 4 < 3 -< i )iftfW 

which provides the basis for the discussion of the scaling properties of Levy flights in an isotropic short range force 
field. 

Deferring details to Appendices B and C we obtain in a precisely analogous manner the renormalization group 
equations for Dx, D 2 , Af , A 2 , A^f , and A 2 , in the case of Levy flights in a general force field with range characterized 
by the indices ql and ar- 



dD x 

~df 

dD 2 

d£ 



= (z- fi)Dx (3.42) 

= (z- 2)D 2 

(Af + Ai)(Di(2 - /*- d) + D 2 (-d)) 



+A 
+A 



[Dx + C 2 ] 2 
(Aj + Al){Dx{2d - 2) + D 2 (2d - 2)) 



[Dx+D 2 Y 
A^d-a L )(Dx+D 2 ) 

+A [Dx+Dtf (3 ' 43) 
dA L A L -I- A L 

+A(2d - 2) i*±m^i (3 . 44) 

[Dx + D 2 \ l 

+A(2d _ 2) w±gm±^i (3 . 45) 

^ = (2,-2-^-4^1^ ,3.46, 

fM = (2 ,_2-„ T )Aj-4^||^ ,3.47) 
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IV. DISCUSSION 



We now turn to a discussion of the renormalization group equations derived i n Se c. III. The Eqs. ( 3.3S| - 3.41 ) 

describe the scaling properties in the isotropic short range case, whereas the Eqs. (3.42-3.47) account for the general 
case of anisotropic force fields with short or long range correlations. We discuss the isotropic short range and long 
range cases in some detail and summarize the results for the anisotropic cases, defe rring deta ils to the Appendix. 

We notice immediately a general feature of the renormalization group equations (3.42-3.47): The requirement that 
the anomalous diffusion coefficient Di, char acter izing the amplitude of the Levy term D\k^, stays constant under 
renormalization, i.e., dDi(£)/d£ = in Eq. ( 3.42 ), immediately implies that the dynamic scale dependent exponent 
z(£) locks onto the scaling index /i, i.e., 



(4.1) 



Consequently, in the case of Levy flights in a weak random force field the long time scaling behavior is entirely 
controlled by the leading anomalous Levy term D^k^ and the dynamic exponent z locks onto /i. In other words, the 
random force field has no influence on the Levy flights. The intrinsic long range superdiffusive behavior, that is the 
occurrence of rare events, enables the walker to escape the inhomogeneous pinning environment and the long time 
behavior is the same as in the pure case. Below we substantiate this claim in more detail when we discuss the isotropic 
short range case. 



A. The isotropic short range case 



In the case of an unconstrained short range force field the renormalization group equations are given by Eqs. ( 3.39- 



3.41 ). They describe how the parameters in the renormalized long wavelength Fokker-Planck equation ( 3.24 ) and 



force correlation function (3.25) change as we differentiably average out the short wavelength degrees of freedom in 
the shell e~ £ < k < 1, characterized by the scale parameter I. 



1. Renormalization group flow and fixed point structure 



Requiring a constant anomalous diffusion coefficient D\ under renormalization, i.e., dD\(£) / dl = 0, Eq. ( 3.39| ) 
implies that the dynamic scale dependent exponent z(£) locks onto /i, z{£) = n, and we obtain the renormalization 
group equations for D 2 and A, 



dD 2 {fi _ 2)D2 + A W : ^+W-2) A 



d£ 
dA 



[Di + D 2 f 
A 2 



e = d c i — d 
d c i = 2/i — 2 . 



(4.2) 

(4.3) 

(4.4) 
(4.5) 

The Eqs. ( [i~2] - |l.3| ) determine the renormalization group flo w in the D 2 — A parameter space. We have introduced 
the parameter e and d c \ is the critical dimension. The Eqs. (4.2-l~3|) determine the renormalization group flow in the 
D 2 — A parameter space. 

Above the critical dimension for e < 0, i.e., d > d c % = 2/z — 2 or [i < l + d/2, Eqs. (4.2-4.3) have the trivial Gaussian 
fixed points D\ = and A* = 0, indicating that (i) the subleading diffusion term, D 2 k 2 , scales to zero compared with 
the leading Levy term and (ii) the quenched disorder, characterized by A, scales to zero and thus is irrelevant. The 
effective long wavelength Fokker-Planck equation takes the form 



(- 



D 1 k fi )P(k,u;)=P (k) 



(4.6) 



and for a particle at the origin at t = 0, i.e., Po(k) = 1, we obtain the scaling expressions in Eqs. (2.5-2.6) with dynamic 
exponent z — /1. Clearly, the physics is characterized by the interplay between the dimension d of configuration space 
and the index /x specifying the tail of the Levy distribution. For /1 < l + d/2 the long range Levy steps predominate 
and control the scaling behavior. 

Below the critical dimension for e > 0, i.e., d < d c \ = 2/i — 2, or 1 + d/2 < /1 < 2, we obtain, solving Eqs. (4.2-1.3) 
for dP 2 /d£ = and dA/d£ = 0, the non-trivial fixed point values for D 2 and A: 
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D* = -Di 



(/i - d)e 



A* = 



(2-/i)(6-2d) + (2-tf)e 

(£>l+£>5) 2 

2A(3 - d) ' 



(4.7) 
(4.8) 



The fixed point D 2 indicates that the subleading diffusive term D 2 k yields a contribution compared to the Levy term 
Dik^. The fixed point value of the diffusion coefficient D\ is negative since the pinning environment created by the 
random force field reduces the ordinary diffusion coefficient from its value D\ — for n < 1 + d/2. The emergence 
of the fixed point A* shows that for d less than the critical dimension d c i = 2/i — 2 the quenched disorder in the 
long wavelength Fokker- Planck equation becomes relevant. In Fig. 11 we have shown the critical dimension d c \ as a 
function of the scaling index fi. For /i = 2we have the Brownian case d c \ = 2; for d < /i < 2 the critical dimension 
d c i depends linearly on fi. Note that d c \ = in the ballistic case for fx = 1. The line d — d c2 — fi delimits the region 
d < d C 2 where naive perturbation theory for 8D2 diver ges as discussed in Sec. III. 

In the Brownian case yj, — > 2 it follows from Eq. (4/7) tha t D\ — > —D\ so that the Levy term Dik* 1 precisely cancels 
the diffusive term D 2 k 2 in the Fokker-Planck equation (3.2); this is consistent with the fact that there is no correction 
to first loop order or more precisely to first order in e = d c \ — din the Brownian case |18|. In Fig. 12 we have plotted 
the position of the fixed point (A*,^) as a function of the scaling index 11 for e = d ci — d > or 1 + d/2 < fx < 2. 



2. Scaling properties 



Introducing the notation P(k, u>, D\, D2, A) = (P(k, u>))f , the scaling properties of (P(k,uj))F are determined in 
the long wave length limit by noting that (P(k, uj)) f ca n equally well be computed from the original Fokker-Pla nck 
equation (3.2) as from the renormalized equation (3.24). From the explicit scaling definitions in Eqs.(3.2C-3.23) we 
thus obtain the homogeneity relation, 



P(k, w, Dx, D 2 , A) = e a ^P(ke e , Dx{t), D 2 (£), A(£)) 



(4.9) 



which determines how (P(k,uj)) p varies as we average out the short wavelength degrees of freedom parametrized by 
the scale parameter £. Similarly, we can derive a homogeneity relation for the wave num ber and frequency dependent 
diffusion c oeffi cient D2(k,u),D±,D2,A), defined according to the Dyson equation (3.4) and k 2 D2(k,u) = — £(fe,u;). 
From Eq. (4.9) we thus obtain 



Da(fc, w, D U D 2 , A) = e 2e - a ^ D 2 (ke e ,Lue a< - l \ D^i), D 2 (£), A{£)) 



(4.10) 



In the vicinity of the fixed point (A*, D*), i.e., for large £, we have, setting from Eqs. (3.30) and (4.1) a{£) — /i£ and 
choosing wave numbers k such that ke e ~ 1, the scaling forms 



P(k, u, D X ,D 2 , A) = k-^Lik/Lu 1 ^, D U D* 2 ,A*) 
D 2 (k,u ) Dx,D 2 ,A) = k^Mik/u 1 ^, A,^,A* 



where L and M are scaling functions. Making use of Eq. ( |3.l| ) we also have 

P{r 1 t,D ll D 2 ,A) - Itl-^GCr/l*! 1 /", A') 



(4.11) 
(4.12) 



(4.13) 



similar to the scaling form in Eq. ( |2.6| ) in the absence of the force field and we infer as in Sec. II a dynamic exponent 
z equal to the Levy scaling index n, i.e., z = fx. 

Above the critical dimension, i.e., for e < or d > d c % =2/j, — 2, we have the Gaussian fixed point (Di,, A*) = (0, 0) 
and we infer from Eq. (4.6) the complex scaling function 



L(x,D u 0,0) 



iDix^ 1 + 1 



(4.14) 
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3. Scaling relations and long time tails 



Using the matching procedure |26||43| , f45| ] we can also derive scaling relations for the w ave numbe r fr equency- 
dependcnt diffusion coefficient D 2 (k,Lu). Expanding the renormalization group equations (U) and Q to first 
order in e about the fixed point (A*,D 2 ) and defining A{£) = A* + SA(l) and D 2 (i) = D* 2 + SD 2 (£) the linearized 
equations take the form dSD 2 /dl — — (2 — n)5D 2 + A((d — n)/D\)8A and dSA/d£ = — |e|<$A with solutions 



SD 2 (£) 



SD" + A 



5A C > 
~~D~i 



-{2-y)t 



SA° 

A e" 

D, 



5A(£) = 6A° e -W 



(4.15) 
(4.16) 



Here SA° = 5A(0) and 8D 2 = 5D2(Q) and Eqs. (4.15-433) hold both above and belo w the critical dimension d c i. 
Since z(€) is locked on to /i and suppressing the dependence on D\ we obtain from Eq. fl4.1C| ) 

D 2 {k, lu, D 2) A) = e i2 ~^ e D 2 (ke e , ue^, D* 2 + SD 2 {£), A* + SA(£)) , (4.17) 

1 and from Appendix B, D 2 (0, 1, D 2 , A) oc A, 



where 5D 2 (£) and SA(£) are given by Eqs. fl4.15|44.16|) . 

In the long wavelength limit k — > and choosing £ such that cue 
we obtain, inserting Eq. ( 1.16| ), the scaling expression 



D 2 (0,uj, D 2 , A) oc oj A* + <5A°^l e l/^ 

Above the critical dimension, i.e., for d > d c \ — 2/i — 2 or e < 0, A* = and we have 

D 2 {0,lu,D 2 ,A) oc 5A Lu {d - d " 2)/d " 2 , 
where we have introduced the other critical dimension 

d c2 = /i . 



(4.18) 



(4.19) 



(4.20) 



In the low frequency limit lo — > the coefficient D 2 vanishes for d > d c2 and diverges for d c i < d < d c2 . This behavior 
is consistent with our remarks in Sec. Ill and clarifies the nature of the divergence. In fact, d c2 — fi plays the role of a 
second critical dimension controlling the behavior of the subleading diffusion coefficient D 2 and supports the heuristic 
argument by Bouchaud et al. in |27[ |. In the Brownian case [i = 2 the two critical dimensions coincide. 
Below the critical dimension, d < d c x or e > 0, A* > and we have the leading behavior 



D 2 (0, oj, Da, A) oc A*uj- {2 - dc2)/d " 2 



(4.21) 



which shows divergent behavior in the low frequency limit. 

From D 2 (t) — J exp(—iLjt)D 2 (uj)duj/2ir we finally obtain the algebraic long time tails for the time dependent 
diffusion coefficient: 



D 2 (0,t,D 2 , A) cx r d/d " 2 for d> d c i, 
D 2 (0,t,D 2 , A) cx r d " ljd " 2 for d< d c i. 

I n a s imilar way we can extract the behavior of D 2 (k, oj, D 2 , A) in the long wavelength limit. Setting to 
(4.17) and choosing £ such that fee 1 ~ 1 we obtain 



D 2 (k,0,D 2 ,A) ocfc-< 2 -^ 



A* + SA°k^ 



(4.22) 
(4.23) 

in Eq. 

(4.24) 



D 2 {k, 0, D 2 , A) oc 6A v k 



i„d-d c 



for d > d c i 



D 2 (k, 0, D 2 , A) oc A*k- {2 - d " 2) for d < d cl 



(4.25) 
(4.26) 



For d > d c2 D 2 vanishes in the long wavelength limit k — > 0; for d < d c2 D 2 is divergent, in accordance with our 
previous discussion. 
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The spatial dependence is inferred from D 2 (r) = J exp(ikr)D 2 (k)d d k/ (2ir) d and we derive the algebraic long range 
fall-off 

D 2 {r, 0, D 2 , A) oc r dc2 ~ 2d for d > d cl (4.27) 
D 2 {r, 0, D 2 , A) oc r^-*-*"- for d < d cl (4.28) 

In Fig. 11 the line d = d c2 , the second critical dimension, delimits the regions for the behavior of D 2 . For d > d c2 D 2 
converges, for d < d c2 D 2 is divergent. 

It is instructive to consider the renormalization group flow in the A — D 2 plane about a fixed point in more detail. 
This discussion is carried out in Appendix D. Another issue, in the analysis of the renormalization group equations 
(3.38)-(3.40) we chose to keep the Levy coefficient D\ fixed under a renormalization group transformation. This 
requirement leads, among other results, to z = /J, which is one of the main conclusions of the present work. Clearly, 
keeping D± fixed is an arbitrary choice and our scaling results cannot depend on this choice. This point is discussed 
in Appendix E. 



B. The isotropic long range case 



We now turn to a discussion of the case of Levy flights in an isotropic long range random force field characterized 
by a fall-off exponent a. The case of Brownian motion in an algebraic long range field has been discussed by several 
authors P, ^0|j2^j49| -^2]| both to first and second loop order. The main conclusion here is that provided the force field 
falls off slowly enough the long range force correlations interfere with the Brownian walk and give rise to anomalous 
diffusion in any dimension. For comparison we have summarized the Brownian case in Ap pendix F. 

Keeping as usual D\ fixed by locking z onto \x we extract from the general equations ( (3.42 - 3.47 ) the appropriate 
renormalization group equations for D 2 and A2 = A, 



dD 2 
dl 
dA 

~dl 



= 0* - 2)^2 
= eA - AA- 



A D x {2d - (i - a) + D 2 (2d - 2 - a) ^ 



A 2 



[Di + D 2 ] 

where we have introduced the expansion parameter 



a c i 



(4.29) 
(4.30) 

(4.31) 



Here a c ± = 2\i — 2 has the same value as the critical dimension defined in Eq. (4.5) for the short range case. In the 
present context a c i, of course, plays the role of a critical fall-off exponent for the long range force correlations. 

The analysis now proceeds precisely as in the short range case. For e < 0, i.e., a > a c i — 2[i — 2 or fx < 1 + a/2, 
we obtain the trivial Gaussian fixed points D 2 = 0, and A* = 0, showing that (i) the subleading diffusion term D 2 k 2 
scales to zero compared with the leading Levy term and (ii) the quenched disorder, characterized by A, scales to 
zero and thus is irrelevant. The effectiv e long w avelength Fokker-Planck equation takes the form in Eq. (4.6) and we 
obtain the scaling expressions in Eqs. ((2.5-2.6) with dynamic exponent z = fi. In contrast to the short range case, 
where the physics is controlled by the interplay between the dimension d of configuration place and the Levy index 
/z, the fall-off exponent a replaces d in the long range case for a < d. For a > a c i or fi < 1 + a/2 the long range Levy 
steps predominate and determine the scaling behavior. 



+ a/2, we obtain the non-trivial fixed point values, D 2 = —D\(2d — /J, — a)e/(4(/x — 
Di,) 2 /AA , indicating that the subleading term D 2 k 2 yields a contribution compared 

-Do, 



For e > 0, i.e., a < a c i or fi > 1 
2) + (2d-2-a)e) and A* = e(D x + 

to the Levy term D\k^ and that the quenched disorder becomes relevant. In the Brownian case fi — > 2 D 2 
i.e., the Levy term Dik^ precisely cancels the diffusive term D 2 k 2 in the Fokker-Planck equation (3.2). 

To leading order in e we have the fixed points (D 2 , A*) = (0, 0) for e < and the fixed points, D 2 = e[(2d — fi — 
a)/4(2— /J,)]Di and A* = eD\/AA for e > 0. Similarly, in the vicinity of either fixed point the linearized renormalization 
group equations d8D 2 jdt— —(2 — fi)SD 2 + A((2d — /i — a)/D\)6A and dSA/d£ — — \e\SA with solutions of the same 
form as in Eqs. ( 4.15| -1.16), i.e., 



Z — /I L)\ 

6A{£) =<5A°e-l £ l f . 



3At-2d-2(5A° _ w 
2- (i D~[ e 



(4.32) 
(4.33) 
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For the distribution function P(k, w, D\ , D 2 , A) and diffus ion c oefficient D2(k, oj, D\ , D 2 , A) we obtain ag ain fo r large 
I and choosing ke £ ~ 1 the scaling forms Eqs. ( 6.67 ) and ( 6.68| ). We also obtain the scaling form in Eq. ( 4.13 ) for P, 
implying that the dynamic exponent z locks onto fi. 

As in the short range case the main conclusion is that in the case of Levy flights in a weak isotropic long range force 
field, the long time scaling behavior is entirely controlled by the leading anomalous Levy term Dik 11 . The random 
force field does not influence the Levy flights. The long range superdiffusive behavior enables the walker to escape 
the inhomogeneous pinning environment and the long time behavior is the same as in the pure case. 

For e < 0, i.e., d > a > a c \, we have the Gaussian fixed point (Dj, A*) = (0,0) and we obtain in particular the 
complex scaling function in Eq. (4.14 ) for P(k, o j). Im plementing the matching procedure as in the short range case, 
we obtain to lowest order in A Eqs. ( 4.18 ) and ( 4.24 ) For e < or d > a > a c \ A* = and we have 



D 2 (0, oj, D 2 , A) oc ,5A w (a - ac2)/Qc 
D 2 (k,0,D 2 ,A) oc(5A fc a - ac2 , 



(4.34) 
(4.35) 



where a c2 = d c2 = /i is the second critical fall-off exponent. In the low frequency limit uj —> or long wavelength limit 
k — > D 2 vanishes for a > a c2 and diverges for a < a c2 . This behavior is in accordance with the naive perturbation 
theory discussed in Sec. Ill and in Appendix B; introducing a force correlation A ~ p a ~ d in the integrand in Eq. 
(3.8) for the correction 5D 2 we obtain a convergent contri butio n for a > a c2 and a diverging 5D 2 for a < a c2 . 

For f > or a < a c i A* ^ and we obtain Eqs. ( 4.21 ) and ( 4.26| ) diverging in the low frequency and long 
wavelength limits, respectively. 

Finally, we obtain for the temporal and spatial behavior of D 2 for e < 



D 2 (0, t, D 2 , A) oc t 
D 2 (r,0,D 2 ,A) cx r 



-a/a c2 



(4.36) 
(4.37) 



and for e > Eqs. ( p~23| ) and ( p~28| ). 

The scaling properties in the isotropic long range case as compared to the short range case are conveniently 
summarized in Fig. 13 and Fig. 14. In Fig. 13 we have plotted the fall-off exponent as a function of fi. The lines 
a = a c2 and a = a c \ delimit three regions I, II, and III. In I the diffusion coefficient D 2 vanishes and the disorder is 
irrelevant, in II the diffusion coefficient D 2 diverges and the disorder is irrelevant, and in III D 2 diverges and disorder 
becomes relevant. We note that Fig. 13 is identical to Fig. 11 with the dimension d replaced by the exponent a. In 
Fig. 14 we have in a plot of the exponent a versus the dimension d contrasted the long range case to the short range 
case. The short and long range cases are delimited by the line d = a. For a > d we have the short range case, for 
a < d the long range case (see Sec. III). In region I delimited by a = a c \ and d = d c \ D 2 diverges and the disorder is 
relevant, in region II delimited by a = a c2 and d = d c2 D 2 diverges and the disorder is irrelevant, and in region I D 2 
vanishes and disorder is irrelevant. In the Brownian case fi — 2 and region II vanishes; the divergence of D 2 coincides 
with the disorder becoming relevant for d = 2 in the short range case and a = 2 in the long range case, see ( |27]]). In 
Appendix B we briefly summarize the relevant aspects of the renormalization group flow in the A — D 2 plane. 

In the case of a general anisotropic force field the fixed point structure becomes more complicated and we encounter 
the same features as in the Brownian case, refs. (J] and [^7|. The analysis nevertheless proceeds as in the isotropic 
case and we therefore defer the discussion to Appendix G (the short range case) and Appendix H (the long range 



V. SUMMARY AND CONCLUSION 



In the present paper we have discussed Levy flights in d dimensions in a variety of quenched force fields with 
fall-off exponent a. The main conclusion of the paper is that the long range characteristics of the random motion, 
characterized by the dynamic exponent z, is essentially not influenced by the quenched force field and the exponent z 
locks onto the step index /i in all cases. In other words, the long range character of the Levy steps enables the walker 
to escape the inhomogeneous pinning environment created by the quenched force field. 

This behavior is entirely different from the case of ordinary Brownian motion where the dynamic exponent z is 
enhanced, corresponding to subdiffusive behavior for d < 2 in the case of short range forces and for a < 2 in the case 
of long range forces. 

Although the dynamic exponent z is unaffected by the quenched environment the phenomenon is still characterized 
by a critical dimension d c \ in the short range case and a critical fall-off exponent a c in the long range case. The 
critical parameters delimit the relevance of the quenched force field. For d > d c \ = 2/i — 2 in the short range case and 
a > d c i = 2/i — 2 in the long range case, the force correlations scale to zero, indicating that the quenched background 
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does not influence the long range walk characteristics. We note that in the Brownian case [i — 2 and, consequently, 
d c i = 2 in the short range case, coinciding with the dimension where Brownian walk becomes transparent, i.e., with 
a finite probability of revisiting a site and thus becoming sensitive to the force correlations. 

For d < d c i or a < a c the strength of the force correlations characterized by A scales to a finite (fixed point) value, 
showing the relevance of the quenched background. In the Brownian case this gives rise to a change of the dynamic 
exponent z, unlike the Levy case where z remains unchanged 

A further aspect of the Levy case is the appearance of a second critical parameter: d c i = H in the short range 
case and a c = \i in the long range case. These parameters delimit the behavior of the subleading diffusive term 
characterized by the ordinary diffusion coefficient D 2 , which, of course, in the Brownian case is the leading term. 
For d > d c i or a > a c the wave number dependent diffusion coefficient D 2 (k) vanishes in the long wavelength limit 
k — > 0; for d < d& or a < a c , correspondingly, D 2 (k) diverges for k — > 0. In the Levy case this behavior is subleading 
and does not affect the leading term characterized by the (unchanged) dynamic exponent z. The above behavior is 
illustrated in Figs. 11,13, and 14. 

In the anisotropic short range and long range cases we encounter as in the Brownian case a more complicated 
fixed point and renormalization group flow structure, depicted in Figs 16 and 17. However, as mentioned above the 
exponent z remains locked onto \x and we find that the behavior of the subleading diffusion coefficient D 2 is the same 
as in the isotropic case. 

Below in Table 1 we have summarized the behavior of D 2 in the various cases. 





D 2 {k) 


D 2 (r) 


D 2 (u>) 


D 2 (t) 


d > 2^ - 2 






ul (.d-fi)/n 


t-d/n 


d<2[i-2 






a; (M-2)/M 


t -(2 M -2)/ M 


a > 2fj> — 2 




r M-d-a 




t -a/u, 


a < 2fi - 2 


k^ 2 


r 2-d-M 


W (M"2)/M 


£-(2^-2)/^ 



TABLE I. Table of the behavior of the wave number, space, frequency and time dependent subleading ordinary diffusion 
coefficient in the case of short range and long range force correlations in the quenched background. 
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VI. APPENDIX 



In the ensuing appendices we include more technical aspects of our calculations. In Appendix A we discuss the 
derivation of a Fokker-Planck equation for Levy flights, in Appendix B we derive self energy, vertex, and force 
corrections to first loop order, and in Appendix C we derive renormalization group equations in the case of Levy 
flights in a quenched force field. In Appendix D we consider the renormalization group flow in the isotropic short and 
long range cases. Appendix E we show that the dynamical exponent does not depend on whether we keep the Levy 
coefficient or diffusion coefficient fixed under scaling. In Appendix F we summarize the isotropic long range case for 
Brownian walks. Finally Appendix G and H are devoted to a discussion of the anisotropic short and long range cases, 
respectively. 



A. The Fokker-Planck equation for Levy flights 



The usual derivation of the Fokker-Planck equation |i|j5J] depends in an essential way on the existence of the 
moments of the distribution. Precisely this assumption breaks down for Levy flights for a step size distribution 
with / < 2 and we must reconsider the derivation of the associated Fokker-Planck equation. In Sec. II we gave 
a heuristic derivation of the Fokker-Planck equation based on the motion in the absence of a force field. Here we 
derive the Fokker-Planck equation generalizing the standard procedure and relaxing the assumption of finite moments. 
Adapting the discussion in ref. [M we define the functional 



/ = J R(r 2 ] 



dP(r 2 \n,t) d 
Ot 2 ' 



(6.1) 



where P(r 2 \ri,t) is the conditional probability distribution and R(r) is a generator function. Using the "chain rule" 
for a Markovian stochastic process, 



P(r 2 |ri,t) = J P{r 2 \r,ti)P{r\ri,t-t x )d d r 
and expanding in Fourier modes, R(r) — J^k elkr Rk, we obtain, replacing dP/dt by (P(t + At) — P(t))/At, 



(6.2) 



I = J2 / d d r 2 e lk ^R k 



dP(r 2 \r u t) 

at 



= J2 d d re lk - r R k 



p ik-(r 2 —r) _ i 

d d r 2 — P(r 2 \r,M) 



P(r\n,t). (6.3) 



In order to evaluate 



ik-{r 2 -r) _ i p ifc.(r(At)-r(0)) _ i 

d d r 2 P(r 2 \r, At) = ( ) 



At v ' J x At 

we invoke the Langevin equation for Levy flights ( |2.11| ), dr/dt = F(r) + rj with incremental solution 

r{At) - r(0) = F(r)At + -qAt . 



(6.4) 



(6.5) 



Inserting Eq. fl6.5|) in Eq.( |6.4[ ), expanding and using the result ( e jfe - I ? A * — 1) = — const. At for a Levy dis trib ution 
(cf. Eqs. ( p^ - ^4| ) we obtain ^ e ^-(F(r)At+r,At _ i)/At) ~ -const.fc^ 1 + ik ■ F(r) which, inserted in in Eq. (|J) and 
requiring it to hold for any variation of R, yields the Fokker-Planck equation for Levy flights (in momentum space), 



dP{k,t) 
di 



= -const. k^P(k, t)-ik-^2 F ( k ~ p) p (p, *) 



(6.6) 
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B. Self energy, vertex, and force corrections to first loop order 



Here we evaluate the self energy, vertex and force corrections to lowest order for Levy flights in a general anisotropic 
force field using the diagrammatic rules and diagrams in Sec. HI. 

The self energy correction is given by the diagrams in Fig. 7. £(fe,a;) splits up in longitudinal and transverse 
contributions, i.e., E(fc,w) = E^fc, u) + Sr(fc,w) . For E^ we have 



£ L (fc,u;) 



-A 



d % ^ {k/2 _ p) K-w-Pnm-P)-im+P) Ga(k/2+p ^ . (6 . 7) 



(2tt 



(k/2-p) 



In order to extract the leading contribution to 0(k 2 ) we first symmetrize Eq. ( |6.7[ ) by replacing p by —p and, further- 
more, expanding G^ 1 {k/2±p, w) and A L (k/2±p) for small k, G " 1 (fc/2±p, w) = Gq x (p, w)±(p • k/2)(nD lP ^- 2 +2D 2 ) 
and A L (fc/2 ±p) = A L (p) ± (p • k/2)(a L - d)p aL ~ d - 2 we obtain to 0{k 2 ) 



1 ^ n A? 

2 



;fc 2 ^ J ^{[(d-2)G -W) + (^i^ + 2£ 2 p^ . (6 

Similarly, we find for Sr(fe, a;) to 0(k 2 ), 



2^ d 



In the static limit w — > and in the isotropic short range case, Al 
particular the expression yielding Eq. (3.8), 



Ay = A, and A L ' T (p) = A, we obtain in 



d d p Di(d - /x)p" + D 2 (d - 2)p 2 



(27r) d 



[£>ip" + L>2J3 2 ] 



212 



(6.10) 



In the context of the momentum shell integration method discussed in Sec. Ill the dilution of degrees of freedom 
requires that the internal momenta fe/2 + p and k/2 — p of the propagator and force contraction, respectively, must 
lie in the momentum shells e~ l < \k/2 +p| < 1 and e~ e < \k/2 — p| < 1. We notice that this assignment is invariant 
under the symmetrization p — > — p performed in order to extract the leading k 2 correction to E. On the shell we thus 
obtain for E(fc,0) in Eq. Qo.lOD for small i, 



E(fe,0) 



1 2 2 S d l D 1 (d-ri+D 2 (d-2) 
2 KA£ ^{2^Yd [D 1+ D 2 ] 2 



(6.11) 



where Sd — 2n d l 2 / (d/2 — 1)! (cf. Eq. (2.1)), leading to Eq. ( 3.32| ) for the correction 5D 2 to the diffusion coefficient. 
In a completely similar manner we obtain in the general case on the shell p = 1, 



E(fe,0) = 

1, 2 S d 1 A| [(d + //)J>i (Af- + Aj) + ( OL - dKA + £ 2 )A^] - A|(2d - 2)(£>! + £> 2 )(Af + A,") g 

~ Ki — — — ™ — £. 



2 (27r) d d 



[Di 



(6-12) 



The vertex corrections are given by the diagrams in Fig. 8. In the isotropic short range case the correction is given 
in Eq. (3.9) and 5X in Eq. ( 3.10D follows directly, i.e., no symmetrization is required here since the correction in Eq. 
(3.9) is already of order k. Focusing on the corrections to Xl and Xt the calculation is also quite simple in the general 
case. We obtain 



SXt 



d d p 



(27T) 



7 A i (p)G (p,0) 



(6.13) 
(6.14) 



and on the shell e 1 < p < 1 
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. a S d 1 A^ + A^ 



(6 - 15) 

Th e force cor rections are given by the diagrams in Fig. 9. In the isotropic short range case the results are given in 
Eqs. ( j3.llH3.12l ). In the general case we obtain the correction to the vertex function T(fc,p, I, uS) defined diagrammat- 
ically in Fig. 6, introducing the tensors (dyadics) A L {k) = [k a k p /k 2 \A L {k) and A T (fc) = [5 al) - k a k p /k 2 }A T (k), 

T(k,p,l,w) = - A|fc • A L {1) -p-\ 2 T k- A T {l)-p 

+ \ 4 L J -^^k ■ A L (n) ■ (p + I - n){k - n) ■ A L (l - n) ■ P G (k - n,uj)G (p + I - n,w) 

f d d n 

+ x t j J2^yl k • A T (n) • (p + Z - n)(k - n) • A T (Z - n) • pG (k - n, cj)G (p + l-n,w) 

/d d n 
- — -rk ■ A L (n) • {p + l - n)(k - n) ■ A T (l - n) ■ pG {k - n,uj)G Q (p + l- n,uj) 
(27r) a 

f d d n 

+ a t a I / 7 — Tl k ' AT ( n ) ■ (p + I - n)(k - n) ■ A L (l - n) ■ pG {k - n, uj)G Q (p + l-n,w) 
J {^) 

/d d n 
^-pfc • A L (n) ■ p(k -n)-A L (l-n)-{p + n)G (k - n, uj)G (p + n, lo) 

/d d n 
^-pfc • A T (n) • p(k - n) ■ A T (l - n) ■ (p + n)G (k - n, oj)G q ( P + n, lo) 

/d d n 
j^yk ■ A L {n) ■ p(k -n)-A T (l~n)-(p + n)G (k - n, lu)G (p + n, w) 

+ A|Ai J-^k-A T (n)-p(k-n)-A L (l-n)-(p + n)Go(k-n,uj)Go(p + n,uj). (6.17) 

In order to identify for example the correction to Aj we choose fe • I = and take the I — > limit. Using that 
k ■ A L (l) p = 0,k- A T (l) p = k-pA T (l),k ■ A T (n) n = 0, and k ■ A L (n) n = k- nA T {n), the expression fl6.17| ) 
reduces considerably and we find the correction 

and on the shell e~ e < p < 1, 



5A± = A£ ( 1 - - ) / j^A 1 (p)A L (p)G 2 (p,0) (6.18) 



sa t_ X 2 S d ( 1\ (Aj + A^)(Af + A£) 

* Ai - Ai Wf V'd) [d 1+ d# / ■ (6 - 19) 

Similarly, we find the correction to A L by choosing k = I, 

5A{ = 5Al , (6.20) 
i.e., the force corrections in the longitudinal and transverse cases are identical. 



C. Renormalization group equations in the case of a general force field 



Here we derive the renormalization group equations for a general force field to first loop order following the procedure 
in Sec. III. Including the correc tions to D -2,\l,Xt, A l and A T we obtain the renormalized Fokker-Planck equation 
and force correlations (cf. Eqs. ( 3.18 - 3~19l )), 
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(-to; + Dxk" + (D 2 + SD 2 )k 2 )P(k, u) = P Q {k) + (A L + SX L ) 



d d p 



(2it) d 



+ (A T + ^A T ) 
(E a (k)E( 3 (p)) F = (A L (k) + 5A L (k)) 
{B a {k)B(\p)) F = (A T (fe) + 5A T (k)) 



{-i)k ■ B(k - p)P(p,uj) 
'k a k ' 



k 2 

§a f3 . 



(2ir) d 5{k + p) 
k a kP 



k 2 



(2ir) d 6{k+p) 



(6.21) 
(6.22) 
(6.23) 



where fe and p after the momentum shell integr ation now range in the interval 1 < k,p < e . Renormalizing wave 
numbers, frequency, and distribution as in Eqs Q3.20H3.22l) and E and B according to E'(k') = exp(-0 L (£))E(k), 
and B(k') = exp(— (3T{£))B(k), where a and (3l,t are given by Eqs. ( 3.30 - 3.31 ) with u re placed by Ul,t we obtain 
a renormalized Fokker-Planc k equation and force correlations of the same form as in Eqs. (3.24-3.25). Furthermore, 
from the expressions in Eqs. ( [2.23 - 2.24 ) for A T (fc) and A L (k) we identify scale dependent parameters and obtain the 
renormalization group equations 

(6.24) 
(6.25) 
(6.26) 
(6.27) 
(6.28) 
(6.29) 
(6.30) 
(6.31) 

1 by choosing 



D[(e) 


= Die" 






D' 2 {£) 


= (D 2 


f SD 2 )e~ 2e 




X' L (£) 


= (Xl- 




-d)£+a+l3 L 


\' T [£) 


= (At - 


f (5A T )e- (1 " 


-d)e+a+/3 T 


AfW 


= (Af 


+ 5A[)e M ~ 




AfW 


= A 2 L e 


U-2f3 L -£(a L 


-d) 


AT(£) 


= (A? 


+ 5Aj)e M ' 


-2/3t 


AfW 


= A T 2 e 


it-1S5 T -t(a T 


-d) 



or in differential form the Eqs. (3.42-3.47). We have here again fixed the vertex constants to Xl,t 
Ul = 1 + d — z — X^SX^/i and ut = 1 + d — z — X f 1 SXt/£- 



D. The renormalization group flow in the isotropic case 



1. The isotropic short range case 



Eliminating I in Eqs. (4.15) and (4.16) and setting [i = d C2 and e = d Cl -dwe obtain 

(2-d C2 )/|, 



5Do 



5D° + A 



5A» 



6A 



SA° 



,-d\ 



A 



6A 



(6.32) 



In Fig. 15 we have depicted the renormalization group flow about the Gaussian fixed point (A* , Z^) = (0,0) in the 
three cases: a) d > d c2 , b) d = d c2 , and c) d c2 > d > d&. In the case a) for d > d C2 , corresponding to region I in Fig. 
11, the first term in Eq. ( |6.32| ) dominates and near the fixed point SD 2 - [5D% +~A5A /D 1 }[5A/5A ]( 2 - d ^/\ d -i- d \, 
i.e., the trajectories approach the fixed point with vertical slope except for the trajectories originating from the 
domain of initial values on the line <5D§ = — (A/ Di)8A° , corresponding to the position of the non-trivial fixed point 
emerging below d c \. In this regime D 2 vanishes as discussed above. In the case b) for d = d C2 , the marginal case, we 
have 5D 2 = (<5Z?2/<5A°)<5A and the trajectories approach the fixed point with constant slope, i.e., the linear scaling 
regime depends on \i and becomes largest for fj, = d, precisely the case where naive perturbation theory or the scaling 
analysis above yield a divergent D 2 . In the case c) for d c \ < d < d c2 , corresponding to region II in Fig. 11, we 
have 6D 2 = —(A/Di)5A. The trajectories now approach the fixed point with a constant slope —A/D\, except for 
trajectories with <5A° = which lie on the line A = . 



2. The isotropic long range case 



Eliminating £ in the renormalization group equations (4.32-4.33) we have 
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5Do 



A- 



6A 
8A° 



(2— d c 2)/|a ci -a 



A^- 2d ~ 2 SA 
2-fj, 



(6.33) 



Focusing on the flow about the Gaussian fixed point (A*, D 2 ) = (0, 0) the discussion in the isotropic short range case 
above applies with d replaced by a. For d > a > /i the trajectories approach the fixed point with vertical slope except 
for the trajectories originating from the domain of initial values on the line 5D 2 = — (A(3/i — 2d — 2) / (2 — /i)-Di) SA°, 
corresponding to the position of the non-trivial fixed point emerging below 2/i — 2. In the marginal case [i = a the 
trajectories approach the fixed point with constant slope, i.e., the linear scaling regime depends on \x and becomes 
largest for /i = a, precisely the case where naive perturbation theory or the scaling analysis above yield a divergent 
D 2 . For a < [i the trajectories approach the fixed point with constant slope — A(3/i — 2d — 2)/(2 — fi)D\, except for 
trajectories with 5A° = which lie on the line A = 0. Fig. 15 also applies in the long range case with d replaced by 
a but note that the slope —A/D\ is replaced by — A(3[j, — 2d — 2)/(2 — fi)D\, 



E. Keeping the ordinary diffusion coefficient constant under renormalization 



We here briefly repeat our renormalization group analysis but now keeping the diffusion coefficient D 2 fixed under 

s calin g. Limiting our discussion to the behavior near the Gaussian fixed point A* = we thus obtain from Eqs. 

( |3. 39 - 3. 41 ) z = 2 and the renormalization group equations dD\jd£ = (2 — [l)D\ and dAl dl = (2 — d)A with solutions 
Di = Di exp((2 — fi)£) and A = A°exp((2 — d)£). From the homogeneity relation (4.9) we have, setting a(£) = 21, 



P(k,u,Di,A) = e 2l P(ke i ,uje 2i ,D°e^-^ e ,A Q e ( - 2 ~' d ^), which forms the basis for our scaling analysis. Choosing 
ke e ~ 1 we deduce the scaling form P(k, u>, D\, A) = k~ 2 P(l,Lu/k 2 , D\k^~ 2 , A°k d ~ 2 ) . Note, however, that this form 
does not imply z — 2 since D\k^~ 2 and A°k d ~ 2 diverge in the long wavelength limit k — > 0. In fact choosing k 
such that A°k d - 2 <C 1 we obtain perturbatively, P(k,uj,D 1 ,A) = k~ 2 (-iuj/k 2 + D\k^~ 2 + D 2 + const. A A; <i " 2 )" 1 or 
P(k,uj,D 1} A) = (— iu + D\& L + D 2 k 2 + const. A°fc d ) -1 similar to our previous results and implying the dynamical 
exponent z = fj>. 



F. The isotropic long range case for Brownian walks 



Here we recover the well-known results in the Brownian case to first loop order [|4|;^7| . The renormalization group 
equations in the isotropic long range Brownian case are extracted from the general equations (3.42-3.47) by setting 



H = 2, Di = 0, D 2 = D, A[ - t = 0, A 



L,T 



A , and qsl.t = a, i.e., 



^R = (z-2)D + A(2d-a-2)^ 



dA 

~dl 



(2z-2- a) A - AA 



D 2 



Keeping D fixed we have z — 2 — A(2d — a — 2)A/D 2 and the equation for A, 

dA , N A A/ A 2 
— = (2-a)A-2A{2d-a)^ . 



(6.34) 
(6.35) 

(6.36) 



The expansion parameter is 2 — a. For 2 — a < we obtain the Gaussian fixed point A* = 0, the disorder is 
irrelevant and z locks onto 2, characterizing ordinary diffusion. For 2 — a > we have the non-trivial fixed point 
A* = (2 — a)D 2 /2A(2d — a) and the dynamic exponent 



z = 2-— — (2 -a) , 



(6.37) 



signalling anomalous diffusion |27|, due to the interference of the quenched long range force correlations with the 
Brownian walk. In the vicinity of either fixed point we have A(^) = const. exp(—|2 — a\£) + A* and we obtain 
scaling relations for the distribution P(k, u>, A) and the wave number and frequency dependent diffusion coefficient 
D(k,tu,A), 



P{k, w, A)= e zl P (ke e , Lue zi , A* + const. e -l 2 - a l 



D{k, u, A)= e^ 2 -^ e D (ke l , uje zl , A* 



const. e 



-a\t^ 



(6.38) 
(6.39) 
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Similar to our discussion in Sec. IV we infer from Eq. ( 6.38| ) the dynamic exponent z and from Eq. (6.39), choosing 
1 ~ 1 and ke l ~ 1, D(0, u, A) oc ^-^1* [A* + co^t.^^l z \ and 0, A) oc k'^^ [A* + const.fcl 2 -* 1 !] 



uJ(t 



. For a > 2 we have A* = and z — 2,, i.e. 

Z?(0,w, A) oc w (q ~ 2)/2 (6.40) 
D(k, 0, A) oc fc a ~ 2 (6.41) 

which vanish in the k — > and a; — > limits, and for the temporal and spatial behavior, 

D(0, t, A) oc t~ a l 2 (6.42) 
D(r,0,A) ocr 2 ~ a - d . (6.43) 

For a < 2 we have A* > and z given by Eq. ( |6.37| ) and 

£>((), w, A) oc A*w~i(^)( 2 - a ) (6.44) 
D(fe,0,A) oc A*/r(^) (2 ~ a) , (6.45) 
where the behavior of D depends also on the dimension of the system d Q^^]. Also, 

D(0, t, A) oc A*i*(fei )( 2 -°)- 1 (6.46) 
D(r, 0, A) oc A*r(^f )(a-o)-d _ ( 6 4? ) 

G. The anisotropic short range case 

Generally the quenched force field must be expected to be anisotropic consisting of a transverse divergence-free 
part B and a longitudinal curl-free part E. The vector character of F is reflected in the transverse and longitudinal 
correlation functions A T and A L , respectively, introduced in Sec. II. The case of Brownian motion in a short range 
anisotropic force field has been discussed to first loop order in refs. |2^|3^]. As in the isotropic case treated in refs. 
^8| , p0| to second loop order, the critical dimension is d c i = 2; below d c \ the long time behavior is controlled by the 
isotropic fixed point A^ = A^ = A*, giving rise to anomalous subdiffusion. However, at intermediate times the 
diffusional character is controlled by a transverse fixed point Ay ^ 0, A^ = 0. 

1. The Brownian case 

In order to clearly illustrate how the Levy case differs from the Brownian case we briefly discuss the renormalization 
group equations in the Brow nian case. Setting D\ = 0, D% = D,A^' T = Al,t, and A^ T = 0, we obtain from the 



general equations (3.42-3.47) 



^ = (2z - d - 2)A T - 2A(3 - d)^^ . (6.49) 

^± = (2z-d- 2)A L - AA^ + 2A(d 1)^1 . (6.50) 
Keeping D fixed, we choose z = 2 — A(2(d — 1)At — dA^/D 2 and we have the equations for At and Ax,, 

= (2 - i)^ - 2A 2 ' J - + j 2 3 - M > A ^ (6.51) 

at D A 

Above the critical dimension d c \ — 2 we obtain the stable Gaussian fixed point (A£, A T ) = (0,0), corresponding 
to normal diffusion for z — 2 and irrelevance of the quenched force field. Below d c \ to 0{2 — d) we obtain the 
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stable isotropic fixed point (A* L , AS,) = [(2 — d)D 2 /2A, (2 — d)D 2 /2A], controlling the long time anomalous behavior 
characterized by z = 2 + 2(2 — d) 2 evaluated to second loop order 0((2 — d) 2 ) (l^JI^] and the unstable anisotropic 
transverse fixed point (AJ, A T ) = (0, (2 — d)D 2 /AA), determining; the cross-over at intermediate times. We shall not 
pursue the Brownian case further here but refer to refs. |2(||35|,^7],fl9] [5^ 1 . 

2. The Levy case 

In the Levy case we p rocee d as in Section IV, setting z = fi in order to fix D±, A^' T — A^y, and A 2 ' T = 0, the 



general equations (3.42-(3.47) imply the renormalization group equations 



-#■ = 0*- 2)D 2 + A (D 1+ D 2 ) 2 (6 ' 53) 

dAL f a o, 2Ai-(d-l)A L A T 

~~d7~ ~ ~ 2A [D 1+ D 2 ] 2 (6 ' 54) 

^=eA T -2A(3-,)^ F , (6,5) 

where e = d c i — d = 2/i — 2 — d. 

Again we observe that irrespective of the vector character of the random force F the dynamic e xponent z locks 



onto the Levy index /i, owing to the long range character of the Levy steps. The equations (6.53-6.55) also have three 
fixed points: 

(D;,A* L ,A* T ) = (0,0,0) (Gaussian) (6.56) 

M>*t>*V = (-*W^y ^2A^y^2l^) ( isotropic ) ^ 

(D* 2 ,A* L ,A* T ) = (| Dl 2 4 _" M 3//) - e ||- ) (longitudinal) (6.58) 

In order to examine the stability of the fixed points we derive renormalization group equations to O(e). Setting 
D 2 (£) = D* 2 +5D(£), A L (l) = A* L +5A L (£), and A T (l) = A^ + SA T (£), we have 

dSD 2 _(^_ 2 + A (d + 2,-4)Al-(2d-2)A* T \ ^ 



di v D i 

, ;t £> 1 (2- M -d) + £>|(2// + d-4) JfA 
D 2 



{2d-2){D 1 -Dl) 

T 
1 



^ ^<5A T (6.59) 



dJAj / 4A* -(d-l)A^A / 2A(d-l )A L 

"dT = I 6 ~ 2A Df ) SAl + { D 2 

dSA? _ ( oA . Q j. A* L \ XA .A T 



e - 2/1(3 - d) — k UA T - 2A(3 _ d)-£ ) JA L . (6.61 



d£ "" V ~ ~'Di) 1 V Df 

In the vicinity of the Gaussian fixed point we have the solution 

^(,) = (^ + A(^ A o + ^ A o)) e - (2 ^ 
A /( 2 -/x-d) XAr0 , 2(d-l) XA0 



„ , , <5AL U + -A '-SA% e tl (6.62) 

D\ \ d — /i d — (J, ) 

SA L (£) = SA° L e ee (6.63) 

tfA R (^) = <5A^e rf (6.64) 

and the fixed point is stable for e < 0, i.e., d > d c i = 2/i — 2. 

Eliminating I we obtain for the flow in the (D 2 , A £ , At) space for e < 
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SD 2 



6D° 



A 
D~i 



A 

2-ii- ,'! 



d — (i 

2(d-l) 



2(d-l] 



6A T 



SA 



L,T 



(2— A*/|e| 



/' 



-«JA, 



SA 7 



(6.65) 



and the discussion in the isotropic case applies with a few modifications. For d > d c2 the trajectories approach the 
fixed point with vertical slope with exception of trajectories lying in the plane 



(6.66) 



In the limiting case d = /x the trajectories approach the fixed point with constant slope, i.e., the linear scaling 
regime depends on \x and becomes largest for \i = d, corresponding to the case where naive perturbation theory yields 
a divergent diff usion coefficient. Finally, for d c \ < ji < d the flow approaches the fixed point tangcntially to the plane 
defined by Eq. ( 6.66| ). The characteristics of the flow are depicted in Fig. 15. 

The scaling analysis of P and D 2 also proceeds as in the isotropic case. Suppressing the dependence on D\ which 
is kept fixed we define P(k, u>, D 2 , Al, At) = (P(k, to)) f and we obtain the homogeneity relations: 



P(k, u, D 2 , A L , A T ) = e^P (k, e e , ue» e , D 2 (l), A L (£),A T (£)) 
D 2 (k, u, D 2 , Al, At) = e^ e D 2 (ke e , we** D 2 (£), A L (t), A T (£)) 



(6.67) 
(6.68) 



For e < 0, i.e., d > d c \,. where the stable Gaussian fixed point controls the scaling behavior, we obtain, of course, 
the same scaling properti es of P as in th e isotropic case, i.e., the dynamic exponent z = fi and P is given by the 
scaling functions in Eqs. (6.67) and (4.14). For the diffusion coefficient D 2 we obtain, correspondingly, 



D 2 (k, w, D 2 , At, At) = eS 2 -^ e D 2 (ke l , cue^, A L e"^ e , A T e~^ (6.69) 

In the long wavelength limit k — > and using the perturbative resu lt D 2 ~ a At + (3 At we have, setting toe^ 1 ~ 1 
and ke^ 1 ~ 1, the same results as in the isotropic case given by Eqs. (4.19) and (4.25) and the corresponding results 
for the temporal and spatial behavior in Eqs. (4.22) and (4.27). 

Below the critical dimension d c \ in the neighborhood of the isotropic fixed point in Eq. (6.57) we obtain from Eqs. 
(3.60) and (6.61) the solutions for 8Al and 8At, 



5A l (£) = - 
5A T {£) = - 



(3-d) 
2(d-2) 

3-d 
2(d- 2) 



5AI 



5A°, - 



?-±SA T 
3-d T 



-W . 



2(d- 2) 

3-d 
2(d- 2) 



1 XA o^„-| £ |^ 



(SA° L - SA° T ) e 
(5A° L ~5A T ) e~ le ^£ 



(6.70) 
(6.71) 



and we conclude that the isotropic fixed pont is sta ble fo r 1 < d < d ci . 
Similarly, near the anisotropic fixed point in Eq. (6.58|) we have 



6A L (t)= JA£ 



±SA° T 



d- 1 



—5A T e~ 



5A T {£) = 5A T e 



(6.72) 
(6.73) 



and we infer that the anisotropic fixed point is stable for d < 1. In Fig. 16 we have depicted the renormalization 
group flow in the (At, A^) plane in the three cases a) d > d c i, b) d c \ > d > 1, and c) d < 1. 

Fin a lly, b elow d c i w e obtain for the diffusion coefficient D 2 the same result as in the isotropic case in Eqs. ( 4. 21 ), 
(p~26] (0|) and (0|). 



H. The anisotropic long range case 



We finally discuss the anisotropic long range case. In Sec. Ill we have given the general renormalization group 
equations in the case of two separate fall-off exponents and ot for the longitudinal and transverse force correlations, 
respectively. For simplicity we here only consider the case of a common fall-off exponent a = ot = cll < d. Keeping 



24 



Di fixed by locking z onto /i and setting A 1 
renormalization group equations 



L.T 



and A 



L.T 



A LiT we extract from Eqs. (0|)-(E4^) the 



dD 2 

di 



= ([i- 2)D 2 + A 



D 1 (2- f i 



Do a 



(D 1 + D 2 y 



-A, 



.D 1 (2d-2) + D 2 (2d-2) . 

-\-A ;— ; — — ; „ fly 



(£>i + D 2 f 



dA L 
di 
dA T 
~df 



eA L - AA 
eA_4A 



(5i + D 2 ) 2 
A L A T 



(D 



(6.74) 
(6.75) 
(6.76) 



where we have introduced the expansion parameter e = 2/i — 2 — a = d c i — a . To first order in e we identify the usual 
Gaussian fixed point (A L , A T ) = (0, 0) and a line of non-trivial fixed points A^ = eDjjAA . We note that here there 
is no isolated non-trivial fixed point unlike in the short range case discussed above p3| . 

To linear order in e we obtain in the vicinity of the respective fixed points, setting A^y = A^ T + SAl.t, 



dSA L 

di 
d5A T 

di 



AA 



SA, 



AA 



— ~ AT 6A T - ^rA* T 5A 



(6.77) 
(6.78) 



SA° L e e 



In the neighborhood of the Gaussian fixed point we have the solutions 5At(£) = SA^e et and 6Al(£) 
showing that the fixed point is stable for e < 0, i.e., for d > a > 2/i — 2 (note that a < d in the long range case). 
Eliminating the scaling parameter i we obtain for the flow in the Al — At plane 5At — (SAj,/SA l )5Al , i.e., the 
trajectories approach the fixed point with constant slope [ |53| ]. Similarly for the scale dependent diffusion coefficient 
D 2 (t) we have 



5D 2 = 



A / 2_j^a SA0 + d(d^l) SA o 



(2-^)1 _ 



2 -»- a 5A« + 2{d - 1] 6A° 



(6.79) 



Eliminating i we obtain for the flow for e < 

'2 — fi — a 



SD 2 



2- pi 



SA\ 



fc^A° 

a — fj, 



SA 



L.T 



$Kt 



(2-/*)/|e| 



A 
~Di 



-SA, 



fi 



2(d-i; 
a — fi 



SAn 



(6.80) 



and the discussion in the anisotropic short range case applies with d replaced by a. 

The sca ling analysis of P and D 2 also proceeds as before. We have the general homogeneity relations in Eqs. 
( |6. 67 -6. 68). For e < 0, i.e., a > 2/i — 2, the Gaussian fixed point controls the scaling behavior. The dynamic exponent 
z = jj, and P is given by Eqs. J6.67fl and ( [4.141 ). For D 2 we obtain 



fceVe^,A L e- |e|£ ,A T e- 



D 2 {k,uj,D 2 ,A L ,A T ) =e- 

and using the perturbative result D 2 oc qAj, + (3At we find, setting we 11 

D 2 (0,uj,D 2 ,A l ,A t ) « a^")/" 



W 



1, and ke e ~ 1, respectively, 



and 



D 2 (k,0,D 2 ,A L ,A T )cxk a -» 

D 2 {0,t,D 2 ,A L ,A T )^t- a ^ 
D 2 { ri 0,D 2 ,A L ,A T )^r^ d - a 

For e > 0, i.e., a < 2/i — 2, we have in the vicinity of the stable line of fixed points SA^t) = 5A° L e~ 



(6.81) 



(6.82) 
(6.83) 



(6.84) 
(6.85) 

ei and 

SA T (£) = <5A° + SA° L (A* T /A* L )(e~ rf - 1) or eliminating i SA T - c5A° = (A T /A* )(SA L - SA° L ) , showing that the 
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trajectories approach the line of fixed points A* L = er,D\/AA with slope A^/A^. In Fig. 17 we have shown the 
renormalization group flow in the At — plane in the cases a) a > d c \ and b) a < d c \ . 
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Figure captions 

Fig. 1. Plot of the Levy distribution for a microscopic step r\. The tail of p(rj) is characterized by the step index /. 
The distribution is cut off at a smallest step distance r]o corresponding to a microscopic length. 

Fig. 2. The scaling index fj, plotted as a function of the step index /. 

Fig. 3. Diagrammatic representation of the Fokker-Planck equation. 

Fig. 4. Diagrammatic representation of the Dyson equation defining the self energy Y,(k,ui) and the renormalized 
propagator G(k,u>). 

Fig. 5. Diagrammatic representation of the renormalized Fokker-Planck equation. 

Fig. 6. Diagrammatic representation of the 4-point vertex function T(k,p, l,co) and the contraction 
/ r(fe, k, I, u)G{k - I, uj)d d l/(2n) d . 

Fig. 7. Diagrammatic representation of the first loop order correction to £(fc,w). 
Fig. 8. Diagrammatic representation of the first loop order correction to A(k,p, u). 
Fig. 9. Diagrammatic representation of the first loop order corrections to T(k,p, 

Fig. 10. Diagrammatic representation of the renormalized Fokker-Planck equation valid in the long wavelength region 
1 < k < e~ e . The slash indicates corrections evaluated on the shell e~ l < k < 1. 

Fig. 11. Plot of the critical dimension d c \ as a function of the scaling index fi. For ^ = 2 we have the Brownian case 
d c i = 2; for 1 < fj, < 2 the dimension d c \ depends linearly on /1. The ballistic case \x = 1 is attained for d c \ = 0. 
The critical dimension d c \ = 2/i — 2 and the line d = /i separate the regions I, II and II, III, respectively. In I the 
diffusion coefficient D 2 vanishes and the disorder is irrelevant, in II D 2 diverges and disorder is irrelevant, and in III 
D 2 diverges and disorder becomes relevant. 

Fig. 12. Plot of the fixed point (A*,D 2 ) as a function of the scaling index fi for 1 + d/2 < [i < 2. 

Fig. 13. Plot of the fall-off exponent as a function of the scaling index /i. For = 2 we have the Brownian case. The 
lines a — ^ and a = 2/i — 2 separate the regions I, II and II, III, respectively. In I the diffusion coefficient D 2 vanishes 
and the disorder is irrelevant, in II D 2 diverges and disorder is irrelevant, and in III D 2 diverges and disorder becomes 
relevant. 

Fig. 14. Plot of a versus d. The line a = d delimits the short range and long range regions. In region III for a > /i 
and d > \x the diffusion coefficient D 2 vanishes and the disorder is irrelevant, in region II for /i > d > 2/i — 2 and 
[i > a > 2\i — 2 D 2 diverges and disorder is irrelevant, and in region III for < d < \i and < a < D 2 diverges 
and disorder becomes relevant. 

Fig. 15. Renormalization group flow in the A — D 2 plane about the Gaussian fixed point (A*,D 2 ) = (0,0). In case 
a) d > fi, case b) d = /i, and case c) d c \ < d < fi. Fig. 16. The anisotropic short range case: Renormalization 
group flow in the At — A^ plane. Case a): d > d c \ = 2/i — 2 and the trajectories flow towards the stable Gaussian 
fixed point (G) with constant slope. Case b): 1 < d < d c \ and the trajectories flow towards the non-trivial isotropic 
fixed point (I). Case c): d < 1 and the trajectories flow towards the non-trivial longitudinal fixed point (A). 

Fig. 17. The anisotropic long range case: Renormalization group flow in the At — A^ plane. Case a): a > 
d c i = 2/i — 2 and the trajectories flow towards the stable Gaussian fixed point (G) with constant slope. Case b): 
a < d c i = 2/i — 2 and the trajectories flow towards the non-trivial line of fixed points (L) with constant slope. 
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